Neutral triplet Collective Mode as a new decay channel in Graphite 
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In an earlier work we predicted the existence of a neutral triplet collective mode in undoped graphene and 
graphite [Phys. Rev. Lett. 89 (2002) 016402], In this work we study a phenomenological Hamiltonian describ- 
ing the interaction of tight-binding electrons on honeycomb lattice with such a dispersive neutral triplet boson. 
Our Hamiltonian is a generalization of the Holstein polaron problem to the case of triplet bosons with non-trivial 
dispersion all over the Brillouin zone. This collective mode constitutes an important excitation branch which 
can contribute to the decay rate of the electronic excitations. The presence of such collective mode, modifies the 
spectral properties of electrons in graphite and undoped graphene. In particular such collective mode, as will be 
shown in this paper, can account for some part of the missing decay rate in a time-domain measurement done 
on graphite. 

PACS numbers: 72.15.Nj72.10.Di 



I. INTRODUCTION 

Recently Novoselov and coworkers have been able to fabri- 
cate graphene, a single atomic layer of graphite 1 . This discov- 
ery has brought graphene to the center of attention of many 
researchers 2 . The fundamental difference of the electronic 
spectrum of graphene with respect to the usual metals is the 
existence of Fermi points around which an effective Dirac the- 
ory describes the electronic states 3 . The suspended graphene 
now can be fabricated in which the effects of impurity and 
substrate is substantially reduced and one can approach the 
ballistic limit of transport with Dirac electrons 4 . 

Starting from a single layer of graphene, and adding further 
layers, one obtains, graphene multi-layers. For few layers the 
even-odd effects due to quantum confinement arise 2 . How- 
ever, as the number of layers exceeds ~ 10, one approaches 
the bulk limit, or graphite. The Dirac part of the energy 
dispersion of graphite is qualitatively similar to graphene 5 . 
The only important difference between the electronic states 
of graphite and graphene is the presence of small pockets up 
to ~ 40 meV, beyond which the Dirac description applies 
to low-energy physics of graphite as well 67 . Ignoring such 
pockets which originate from the weak interlayer coupling, 
the electronic structure of bulk graphite can be approximately 
described by a tight binding model on a 2D honeycomb lat- 
tice. In our approach both highly oriented pyrolitic graphite 
(HOPG) as well as undoped graphene are treated within this 
model. The calculations of this paper is aimed to explain the 
life-time anomaly in HOPG, but applies to undoped graphene 
as well. 

The presence of Dirac points makes the nature of particle- 
hole excitations in graphene, drastically different from sys- 
tems possessing extended Fermi surface. Due to such cone 
like spectrum, there will be a region below the particle-hole 
continuum, where no particle-hole pairs can exist. Such a 
"window" does not exist in usual metals 8 . A simple ran- 
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dom phase approximation (RPA) analysis shows that presence 
of such window below the particle-hole continuum, provides 
a unique opportunity for existence of a triplet bound state 
of electron-hole excitation 8 . An intuitive way to understand 
such a triplet electron-hole bound state is to view the semi- 
metallic graphene from semiconducting side. From this point 
of view, such collective excitation can be regarded as analogue 
of triplet excitons 9 . 

In this work we focus on the life time experiment done on 
HOPG sample which corresponds to undoped graphite. The 
time resolved photoemission spectroscopy (TRPES) done by 
Moos and coworkers 10 on HOPG, was employed to measure 
the decay rate of quasi particles in graphite. There are two 
salient features of the TRPES experiment reported by Moos 
et. al. 10 which for convenience has been included in Fig. 1 : (i) 
The plateau in the energy range 1 — 2 eV is already a marked 
deviation from Fermi liquid prediction which was qualita- 
tively explained in Refs. ' , in terms of the peculiar from of 
the graphite dispersion near the saddle point. Such a plateau 
has been reported in the carrier life time of doped graphene 
in angular resolved photoemission spectroscopy (ARPES) ex- 
periments as well 12 which can be understood in terms of a 
similar GoW type of treatment 13 , (ii) The second important 
observation of the above TRPES experiment was that the de- 
cay rate of excitations in the whole range of energies over 
which the measurement was performed, was larger than the 
ab-initio calculation of Ref. 1 1 . This clearly means that there 
should be another decay channel for quasi particles, especially 
in the energy range 1 — 2 eV. In the whole measurement range 
the experimentally observed decay rate is almost a factor of 2 
larger than the GW calculation. 

Obviously the phonons cease to exist beyond 0.2 eV, and 
hence can not be responsible for the missing decay channel 
in the energies reported in Ref 10 . Moreover, both in HOPG 
and undoped graphene, there are no plasmons whatsoever 14 . 
Therefore we believe this lifetime experiment already point 
to the existence of an unnoticed bosonic branch of neutral 
excitations 8,9 . There are also other evidences based on the 
Fermi velocity renormalization measurements: If one appeals 
to electron-phonon coupling to explain the experimentally ob- 
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served reduction in Fermi velocity vp with respect to band 
structure prediction, one has to use an electron-phonon cou- 
pling which is almost ~ 5 times larger than the density func- 
tional theory estimates 17 . Therefore it seems that the phonons 
are not enough to account for about 20% Fermi velocity renor- 
malization 17 . The second experimental hint for the existence 
of such a bosonic mode, is the remarkable observation of the 
Bose metal-insulator transition tuned by magnetic field 18 . 

Based on the above evidences, in this paper we employ a 
triplet bosonic mode predicted in Ref. 8 with a gapless disper- 
sion up to uj max ~ 2.1 eV. Our model is a natural general- 
ization of the polaron problem, with spin-flip processes in- 
cluded. We generalize the momentum average (MA) approx- 
imation developed in the context of the polaron problem by 
Berciu 19 to take into account the spin-flip vertices as well as 
the nontrivial dispersion in the spectrum of bosons. First we 
introduce our model and the MA method. Next we apply the 
MA approximation to discuss the coupling of a triplet boson 
to electronic states of graphene quasi particles. The details of 
generalization of MA approximation to spin-flip processes is 
discussed in the appendix. 

n. MODEL AND METHOD 

We start with the Hamiltonian (1) describing the tight- 
binding electrons on honeycomb lattice (first term), along with 
dispersive triplet bosons (second term) and the interaction be- 
tween electrons and bosons (third term): 

£,as=T,l q,m=±l 

E ( s L + s -^'H- ? , a % a n (i) 

where = ±tyl 1 + cos (V3ky/2) cos(fc T /2) + 4cos 2 (fc x /2) 
is spectrum of fermions for (conduction/valance) band and lo^ 
describes the dispersion of spin-1 bosons 8 . Here ct (eg Q ) is 
creation (annihilation) operator for fermions with momentum 
k and spin a =f , J, in either of the valence or conduction 
bands, while St m ,S^ m are ladder operator for spin-1 
bosons with momentum q, and magnetic quantum numbers 
to = ±1,0. In this Hamiltonian, g is the coupling strength 
and describes how strongly the exchange of triplet excitons 
takes place among the electrons. Estimates of a similar 
coupling in doped solid Cgo suggests g ~ 0.3 for those 
systems 20 . Presence of such term, favours singlet pairing 
under suitable conditions 20 ' 21 . 

The interaction term of the Hamiltonian (1) describes both 
spin flip (m = ±1) as well as non spin flip (to = 0) processes. 
Since non spin flip processes can exist in presence of spin- 
bosons as well, to isolate the contribution of the spin flip 
processes, we focus on m = ±1 terms only. In this sector, 
requiring the Hamiltonian (1) to be Hermitian, gives rise to 
the following restrictions on possible values of m, ml ', a, a': 

m = 1(-1) - to' = -1(1) - a =1 (T) - a' =| (|). 



We use momentum average (MA) approximation to calcu- 
late the Green's function and self-energy of system 19 22 which 
yields various physical quantities such as decay rate. Com- 
parison of MA and its descendants (e.g. MA(1), MA(2), etc) 
with other methods demonstrated that this method is accurate 
for the entire spectrum (both low and high energy) and for all 
coupling strengths and in all dimensions 22 . This approxima- 
tion was also used successfully for analysis of the effects of 
ripples on graphene sheet 23 . In the following, we use MA(1) 
approximation, the details of which for the case of dispersive 
mode with spin-1 are derived in the appendix. 

The single electron Green's function can be written as: 

G a>0 {k,T) = -i^rXOlc^e^cyO), (2) 

where a, j3 are spin indices, and |0) is vacuum state. In the 
absence of bosons the free propagator is 

G (k,u>) = (3) 

U — + If] 

To take into account the coupling to triplet bosons, we use the 
equation of motion for G a< p(k, r) to obtain (see appendix), 

G a ,p{k,u>) = G {k,u>)[6 a ^ +g^2F°~P(k, ^,mi;u)], 

91, mi 

(4) 

where 

uj — H + ir] K q - p 

(5) 

Here, F\ is the amplitude for the process in which the ini- 
tial state contains a fermion and a boson, and the final states 
contains only a fermion with opposite spin. Hence, physi- 
cally it corresponds to the amplitude of annihilating one triplet 
(Am = ±1) boson. Applying again the equation of motion 
to Fi generates hierarchy of equations containing amplitudes 
with multi-boson states: 

F"'~ (k,qi,mi,u>)= (6) 
G (k - qi,u - u>$ 1 )[g 2 +^Fg'P(k,q\,q2,mi,m 2 -,u>)]. 

92, "12 

Although each internal vertex may contain spin flip scatter- 
ings, since the Hamiltonian (1) preserves the spin, the incom- 
ing and outgoing fermions must have the same spin. Hence 
the Green's function (2) must be diagonal with respect to the 
spin indices. The rigorous proof of this is given in the ap- 
pendix. Also, by Dyson equation, the self-energy is also diag- 
onal with respect to spin indices: 

G a , a (k, u;) = [Lu-ej:~ £ Q ' Q (fc, u) + ir,]" 1 ■ (7) 

The self-energy S a,a (fc,w) in MA(1) approximation is 
given by (see appendix), 



3 



5 2 EjE, gi G o(>-<fi,w-w ?1 -5 2 Ai(w)) 



1 - .9 2 EjE.ii G (fc - fi, w - w fi - 5 2 AiH)(A 2 H - ' 



(8) 



where A\ , A2 are functions of u>, defined in the appendix. The 
self-energy contains all interaction effects, and can be used to 
calculate spectral weights, decay rates, etc. in a straightfor- 
ward way 24 . 



III. RESULTS 

Now we are in position to derive decay rate or lifetime 
of quasi particles (QP) of HOPG/graphene in presence of 
spin-1 bosonic collective mode. There are some other de- 
cay mechanisms such as electron-hole 11 , electron-phonon 16 , 
and electron-plasmon scatterings 14 . In doped graphene, all 
the above mechanisms might contribute to the renormaliza- 
tion of QP properties. However, in HOPG graphite and un- 
doped graphene there are no plasmons to couple to electronic 
degrees of freedom. 




A. The decay rate 

The imaginary part of self-energy related to life-time and 
decay rate of QP, 

- oc -Im[TrS Q ' Q (fc,w)]. (9) 

T 

We have numerically evaluated the integrals necessary to 
get the self-energy (8). In Fig. 1, we have plotted the de- 
cay rate measured in TRPES experiment of Ref. 10 (trian- 
gles) along with the electron-hole decay mechanisms captured 
within GW approximation (open circles) 11 . As can be seen 
in this figure, the decay into incoherent electron-hole pairs 
within GW approximation can only account for half of the ex- 
perimentally reported QP decay rate. In this figure, we plot the 
total decay rate in presence of the spin flip scatterings from a 
tripled bosonic mode for the coupling values g = 0.25 (filled 
circles) and g = 0.3 (dashed line). The triplet bosonic collec- 
tive mode has a wide dispersion between zero to u) max ~ 2.1 
eV. 

As can be seen, a dispersive bosonic collective mode can 
account for the missing decay rate in HOPG graphite. The 
same result applies to undoped graphene as well. The fact 
that GW approximation falls a factor of two behind the ex- 
perimentally measured decay rate, indicates that in addition 
to incoherent electron-hole decay processes, there should be 
another decay channel provided by a coherent bound state 
of electron-hole pairs, which is what our phenomenological 
Hamiltonian (1) describes. A simple RPA analysis showed 
that such a bound state can occur in triplet channel 8-9 . 



FIG. 1: (Color online) Quasi particle decay rate in HOPG graphite. 
The triangles indicate TRPES measurements in Ref. 10 , the open cir- 
cles are ab-initio GW calculation of Ref." . Filled circles and dashed 
lines show the total decay rate in presence of a new decay mech- 
anism caused by triplet bosons for the electron-boson couplings 
g — 0.25, 0.3. In this calculation we have taken w max = 2.1 eV. 
See the text for explanation. 

B. Dependence on cj ma x 

The dispersion of triplet bosonic mode is over a wide en- 
ergy range from zero to ui max ~ 2.1 eV. The shape of dis- 
persion and the value of w max in the original work or Ref. 8,9 
is essentially controlled by the short range part of the interac- 
tion (Hubbard U). It was also shown that the long-range part 
of the Coulomb interaction does not play a crucial role in the 
dispersion of the spin-1 collective mode 9 . In the present cal- 
culation, we have fitted the dispersion relation obtained from 
the RPA analysis of Refs. 8 - 9 with ~ 10 cosine harmonics, over 
the whole Brillouin zone. 

In Fig. 2 we explore the dependence of decay rates on the 
dispersion bandwidth (w max ). Left panel shows the imaginary 
part of the self-energy for various values of the electron-boson 
coupling g, and for u> max = 1.4, while the left panel shows the 
same result for co max = 2.1. As can be seen in both panels, by 
increasing the coupling strength g, the decay rate at a given 
energy scale increases. Comparison of the left and right pan- 
els for the same values of g shows that smaller width of dis- 
persion (w max ), the bosonic mode leads to stronger spin flip 
scattering. The limit u> max — > can be thought of an Einstein 
like phonon mode which was studied within MA(1) in Ref 23 . 
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FIG. 2: (Color online) Imaginary part of self-energy vs energy (eV) 
for uj(q)max = l.4eV,uj(q)max = 2.1eV. 



Smaller w max in our phenomenological Hamiltonian (1) cor- 
responds to larger U in the Hamiltonian of the original elec- 
trons in Ref. 9 . Hence the observation of Fig. 2 can be justified 
as follows: In terms of the Hubbard type Hamiltonian of Ref. 8 , 
larger U naturally leads to stronger decay rates. 
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FIG. 3: The spectral weight along high-symmetry cut in the BZ for 
Umax = 2. lei/ and (a) g = 0.3 for different energy from -0.6 to 2.6 
eV, (b) g = 0.6, (c) g = 0.8, (d) g = 1.0 for different energy from 
0.2 to 2.6 eV 



C. Spectral function 

Once we calculate the self-energy £ Q,Q (u>) at any approx- 
imation, we are able to immediately calculate the spectral 
weight A(k, u>) = — ^Im[TrG a>Q: (fc, u)]. We have plotted 
the spectral weight along high-symmetry cut r — K — M of the 
Brillouin zone in Fig. 3. We have plotted the spectral weight 
for different energies. Panels (a)-(d) correspond to different 
values of g as indicated in the figure caption. 

The first point to notice in all panels is that the cone like 
dispersion of the Dirac electrons remains quite robust against 
increase in the electron-boson coupling g. To see this more 
transparently, in panel (a) we have plotted some negative en- 
ergy spectral functions as well. As can be seen in panel (d), 
large values of coupling g ~ 1 lead to a remarkable broad- 
ening in the quasi particle peaks. Direct comparison with 
ARPES experiments on graphene indicates that g can not be 
as large as g ~ 1. 

Negative energy plots of panel (a) indicates that there is an 
asymmetry between the positive energy states and the negative 
energy states. This is natural, as the collective mode is an 
excitation and does not carry negative energies. 



undoped graphene . Employing the momentum average self- 
energy we showed that such a coherent particle-hope bound 
state in triplet channel can account for substantial part of the 
missing decay rate in TRPES experiment of Ref. 10 in HOPG. 
Another supporting evidence for existence of such a spin-1 
collective mode which is a natural generalization of triplet 
excitations of ordinary semiconductors to the case of semi 
metallic HOPG comes from the downward renormalization of 
vf 11 ■ Apparently phonons fail to account for such renormal- 
ization. Moreover, the remarkable observation of Bose metal- 
insulator transition tuned by magnetic field in HOPG 18 , might 
indicate that there such spin excitation branch can have inter- 
esting consequences for the behavior of HOPG and graphene 
in magnetic fields. 
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and for odd N: 



(a) 



<b) 



(c) 



(d> 



FIG. 4: First and second order Feynman diagrams for the scattering 
vertex from a spin-1 collective mode. At each vertex the spin of 
electron is flipped. Therefore, the incoming and outgoing spins end 
up to be identical as in diagrams (a), (c). However, if we insist to 
have the spin of outgoing state to be opposite to that of incoming 
state, some vertices (denoted by <g>) will not do flip the spin; that is 
they will not correspond to scattering from a spin-1 collective mode, 
as in diagrams (b), (d). 



E^ J? v..,^ + i^+i(^9i,92,...,+i,-i,...;w) = 

B(k,q 1 ,q 2 ,...,+l,-l,...;Lj)G r , L (k,iv)=0. (A8) 
The A, B coefficients for various orders can be seen by in- 
spection to be non zero. This proves that the spin off-diagonal 
components of the Green's function are zero. This can be seen 
intuitively in Fig. 4. 

To proceed further, we define a modified form of the 
bosonic Green's function in MA(1) approximation as: 

fn~ a (k,qi, -,q n , ■■■ + 1,-1, -M = 

g n Fn' i ~ 1 )a (fc,gi,...,g;,... + 1,-1, ..-M 
G at0l (k,ui) 

By inserting this in Eq.(4) one finds, 



(A9) 



APPENDIX A: GENERALIZATION OF MA(1) FOR 
SPIN-FLIP HAMILTONIANS 

We start by writing Eqn. (4) with explicit spin indices. The 
matrix elements of Green's function become, 

G T , T (fc^)=Go(fc,^)[l+ 5 ^-^ 1 Ta (fc,g 1 ,+l;^Al) 
G ia (fc,w)=G (fc,w)[l+ 5 Efe, fl ^i'' T (^5i 5 -l;^)](A2) 
G u (fc,w)= G (M)bEk, 9l *i' T (*>?i,-l;w)] (A3) 
G LT (fc,w) = G (*,w)[sEfc lS i Ji U (*»&.+l;w)]. (A4) 



G ata (k,Uj) = 

Go(fc,w)[l+ 5 ^/ 1 Q '- Q (fc, (7 l,l;c I ;)G ct , C[ (fc,a;)].(A10) 

k,qi 

where the matrix form of the Green's function is: 

GM -(«*•#-> 0u ° t , B) ). (A") 
Dyson's equation, 



As can be intuitively seen in Fig 4, the non diagonal ele- 
ment of Green's function should be zero. To see this more 
systematically, one writes the one bosons Green's function as, 



^if T (fc,gi,-l: 

k,qi 



gG (k,oj) 



(A5) 



Repeating the equation of motion we obtain the two boson 
amplitude: 



J2F^' l (k,q 1 ,q 2 ,-l,+l;ij) = 



(A6) 



gG (k,uj) 



gG Q (k -q,uj~ u)(q)) 



_gG (k,uj) 

Finally for order N + 1, we obtain for even A^: 

52n&&,...,fr+i F n+i$, 91,92, +1, -1, ...;w) = 
A(k,q\,q 2 ,..., +1,-1,.. .;u)G hl (k,u) = 0, (A7) 



G a . a (k, u) = [u-e % - E a ' a (k, oj) + in}' 1 , (A12) 
give the spin diagonal self-energy as, 

S Q ' Q (uj) = J2 fi'~ a (k, (A13) 

k.qx 

So the MA(1) self-energy is diagonal and can be casted into 
the final form given by Eq. (8), where 



A 2 (u) 



2go,3(^)go,2(^) ' 



2 ffo,2(w) 



3go,3(^)go,2(^) ' 



1 - 



(A14) 



1-. 



fe,«l,...,fn 
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